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It is shown that the kinematic interaction caused by the quasi-Fermi character of commutation 
relations for operators of the atomic representation can induce pseudogap behavior of the spectral 
characteristics of an ensemble of Hubbard fermions. Mathematically, the presence of the kinematic 
interaction manifests itself in modification of the faithful representation of a single-particle Green's 
function of Hubbard fermions D(k, ioJn), which involves, apart from self-energy operator ^^(k, iuJn), 
strength operator P(k, icjn). It is important that the strength operator enters both the numerator 
and the denominator of the exact expression for D(k,iujn)' The kinematic interaction, therefore, 
not only renormalizes the spectrum of elementary excitations but significantly affects their spectral 
weight. It results in strong modulation of spectral intensity A(k,u}) occurring on a Fermi contour. 
Calculations of the spectral properties for the t — J model in the one-loop approximation yield good 
quantitative agreement with the ARPES data obtained on cuprate superconductors. 

PACS numbers: 71.10.Fd, 71.18.+y 71.27.+a, 74.40.-n, 74.72.Kf 



I. INTRODUCTION 

The normal phase of cuprate superconductors is char- 
acterized by a number of intriguing properties that can- 
not be described within the traditional Fermi liquid the- 
ory. These are, first of all, the pseudogap behavior in 
an undoped region^^. It still has been unclear whether 
this behavior is explained basing on the modified Fermi 
liquid concept or requires building the ground state of a 
new type"^^^. In the phase diagram of cuprates, the region 
of parameters corresponding to the pseudogap state bor- 
ders on the region of the superconducting state. This fact 
is often interpreted as resemblance of the formation of a 
pseudogap and of Cooper instability. Since in many stud- 
ies the nature of coupling in cuprate superconductors is 
attributed to spin fluctuation processes, these processes 
are considered to be the fundamental cause of the for- 
mation of the pseudogap state. In recent years, study of 
the interrelation between electron and spin subsystems 
has become especially important for understanding the 
formation of the pseudogap state and establishing the ef- 
fect of pseudogap peculiarities of the spectral properties 
on Cooper instability. This has been a subject of numer- 
ous experimental and theoretical works on physics of the 
strongly correlated systems. 

The phenomenological concept of the formation of the 
pseudogap state proposed in study^ uses the idea of the 
formation of a quantum spin liquid following the scenario 
of the resonance valence bond (RVB) method developed 
by Anderson^. This approach allowed reproducing the 
formation of a pseudogap, growing with decreasing dop- 
ing level, near the antiferromagnetic Brillouin zone. This 
process is accompanied by rearrangement of the large 
Fermi surface experimentally observed in the optimal 
doping region to the small Fermi or Luttinger pockets 
arising in the undoped region. This model was success- 



fully used in the description of the effect of the pseudogap 
on certain properties observed in undoped cuprates, such 
as electronic specific heat^, London penetration depth^, 
superconducting gap^^, electron density of states^^, elec- 
trical and heat conductioni^, and anomalies of optical 
conduction^^. The model was also employed to interpret 
the photoemission spectroscopic data (ARPES) 

The pseudogap formation was also studied using nu- 
merical calculations on the basis of exact diagonaliza- 
tion^^, the quantum Monte Carlo method^^^^^, and the 
dynamic mean field theory (DMFT)i^ii^. The results 
obtained were in satisfactory agreement with the exper- 
imental data, specifically, the presence of a large Fermi 
surface at optimal doping^^ and modulation of spectral 
intensity and reduction of the density of states at the 
Fermi level at weak dopingi^ii^i^. 

Analysis of the spectral properties of cuprates is of- 
ten based on the Hubbard model^^ and its low-energy 
version, i. e., the t — J model. Study of the Hubbard 
model by the renormalization group method revealed de- 
viation of its spectral properties from those described 
within the theory of an ordinary Fermi liquid and a no- 
ticeable decrease in the Fermi surface area^^^^^. The 
features of the pseudogap state directly related to an- 
tiferromagnetic spin fluctuations (SFs) were considered 
in studies^^'^^ within the phenomenological spin-fermion 
model^- and by phenomenological investigation of the 
Hubbard model^^. In study^^, using the method of equa- 
tions of motion^^ and Mori's technique of projection op- 
erators^^, the spectral functions and Fermi surface were 
investigated in the framework of the t — J model. It 
was shown that long-wavelength SFs govern the low- 
frequency behavior of a system, leading to truncation 
of a large Fermi surface, modulation of spectral inten- 
sity, and a decrease in the density of states at the Fermi 
level in the weak doping region. A microscopic theory 
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for the electron spectrum of the Cu02 plane within the 
Hubbard model was proposed in study^^. In this study 
the Dyson equation for the single-electron Greens func- 
tion in terms of the Hubbard operators was derived and 
solved self-consistently for the self-energy evaluated in 
the noncrossing approximation. Electron scattering on 
spin fluctuations induced by the kinematic interaction 
was described by a dynamical spin susceptibility with a 
continuous spectrum. At low doping, an arc-type Fermi 
surface and a pseudogap in the spectral function close to 
the Brillouin zone boundary was observed. 

In a number of studies on microscopic investigation of 
the spin-fluctuation nature of the pseudogap state, the 
interaction between electrons and a spin density wave 
was used as a mechanism of the spin-electron correla- 
tion. The existence of such a wave was considered to 
be an a priori specified property and its origin was not 
discussed. Meanwhile, by now the existence of the spin 
density wave in cuprate superconductors has not been 
experimentally confirmed. In view of this, a reasonable 
question arises concerning possible implementation of the 
pseudogap phase at the interacting electron and spin de- 
grees of freedom but with no use of the hypothesis of the 
spin density wave existing in the system. To answer this 
question, one should take into consideration the impor- 
tant feature of strongly correlated systems, which include 
high-temperature superconductors. As is known, in the 
regime of strong correlations, the adequate description of 
electron systems that takes into account Hubbard corre- 
lations is based on the atomic representation^^. In this 
case, the operators employed in the theory do not satisfy 
the commutation relations characteristic of the Fermi op- 
erators, since commutation of two basis operators of the 
atomic representation results in a basis operator and not 
in a number 

[Xf , = 5f^{5,rXf^ ± S,,X}^). (1) 

Physically, this feature of the commutation (kinematic) 
relations between basis operators manifests itself (for in- 
stance, in derivation of equations of motions or calcu- 
lation of a scattering amplitude) as an additional inter- 
action arising in a system. This interaction, with re- 
gard to its nature, is named kinematic. The occurrence 
of the kinematic interaction in Heisenberg ferromagnets 
was mentioned by Dyson [33] . This interaction occurs in 
them due to the noncommutative character of the SU (2) 
algebra of spin operators. Taking into account the kine- 
matic interaction, Dyson performed the correct calcula- 
tion of a two-magnon scattering amplitude and obtained 
valid temperature renormalizations for both the elemen- 
tary excitation spectrum and the thermodynamic char- 
acteristics [33]. 

In cuprate superconductors belonging to the Hubbard 
strongly correlated systems, the kinematic interaction 
manifests wider. Apart from renormalizing the proper- 
ties of the normal phase, this interaction can be a mech- 
anism of Cooper instability^. The kinematic interaction 
originates from the fact that, in the regime of strong cor- 



relations, the Hubbard model is adequately described on 
the basis of the atomic representation. In this case, the 
Hubbard operators are basis. For them, commutation 
relations are more complex than those for spin opera- 
tors, since they include both quasi-spin and quasi-Bose 
operators. Commutation of two quasi- Fermi operators 
results in a quasi-Bose operator expressed via the quasi- 
spin operators and operators reflecting charge fluctua- 
tions. Thus, the kinematic interaction in the systems of 
interest couples Hubbard fermions with charge and spin 
fluctuations. 

The dynamic and kinematic interactions of Hubbard 
fermions lead to renormalization of their energy spec- 
trum. Therefore, the faithful representation for the 
single-particle Green's function of Hubbard fermions 
DCk^iUn) involves, apart form self-energy operator 
^^(k, icjn) caused by the dynamic interaction of Hub- 
bard fermions, strength operator P(k, icj^) arising due 
to the kinematic interaction of these fermions. It is im- 
portant that P(k, icJn) enters both the numerator and 
the denominator of the faithful representation for the 
distinguished Green's function. While the occurrence of 
the strength operator in the denominator of I^(k, zcj^) 
directly affects renormalization of the elementary excita- 
tion spectrum, the occurrence of the strength operator in 
the numerator of I)(k, zcj^) determines renormalization 
of spectral intensity. This is of fundamental importance 
for investigation of the spectral characteristics of strongly 
correlated systems. 

The above-mentioned coupling of the quasi-Fermi op- 
erators of the atomic representation with the quasi-Bose 
operators implies that, physically, the kinematic inter- 
action reflects the presence of the interaction between 
charge and spin degrees of freedom and between Fermi 
and Bose excitations in a strongly correlated electron 
system. Therefore, the calculation of contributions of 
these interactions to renormalizations of energies of the 
elementary excitations and their spectral intensities is re- 
duced to the calculation of the self-energy 5]L(k, iun) and 
strength P(k, iun) operators. This is the specific way of 
theoretical study of the pseudogap phase in a spatially 
homogeneous case. The present study is aimed at solving 
this problem within the t — J model on the basis of the 
diagram technique for Hubbard operators^^-^. The key 
point of the developed theory is that it uses the faith- 
ful representation for the Matsubara Green's function 
I)(k, zcJn) via the self-energy ^^(k, zcj^) and strength 
P(k, zcjn) operators. 

In Section [III using the modified Dyson equation for a 
strongly correlated system, the correlation between the 
spectral intensity and the strength and self-energy oper- 
ators for the t — J model is established. Then, using the 
diagram technique for Hubbard operators, contributions 
to the strength and self-energy operators are calculated 
in the one-loop approximation and the integral equation 
for the correction to the strength operator is written. In 
Section Hill the choice of magnetic susceptibility entering 
the integral equation kernels is discussed. In Sectior flVl 
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the results of the calculation of the Fermi excitation spec- 
trum, spectral intensity, Fermi surface, and density of 
electron states demonstrating the pseudogap behavior of 
the system are reported. Section |Vl presents the data for 
the limiting case of strong electron correlations, which 
allows obtaining relatively simple analytical expressions 
for the energy of Fermi excitations and spectral intensity. 
The kinematic mechanism of the pseudogap state forma- 
tion is demonstrated in the microscopic scale. The final 
section contains discussion of the results. 



II. CORRELATION BETWEEN THE 
SPECTRAL INTENSITY AND THE STRENGTH 
AND SELF-ENERGY OPERATORS 



The Hamiltonian of the t 
representation is 



J model in the atomic 



H -- 
Ht 



(2) 
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where X^^ = \f^p){f^q\ are the Hubbard operators^^ de- 
scribing the transition of an ion in the f-th site from the 
one-site state \f,q) to the state e is the energy of 

one-electron one-ion state, /i is the chemical potential of 
the system, a = ±1/2 (a = —a) is the spin moment 
projection, tfm is the integral of electron hopping from 
the m-th to f-th site, Jfm = 2t'j^/U is the exchange 
integral, and U is the Hubbard repulsion parameter. 

We calculate spectral intensity A(k, uj) with the use of 
the diagram technique for Hubbard operator o^^i^^ , intro- 
ducing the Matsubara Green's function 



OaMf.r-J^r') = -{TrX'f^r)Xf{r')) = (3) 

]<.UJn 

Here is the operator of Matsubara time ordering. In 
Expression ([3]), the Hubbard operators are taken in the 
Heisenberg representation with Matsubara time r 



vO' 



'{t) = exp(ri^)X^^exp(-ri^), < r < 



1 



(4) 



where T and H are the temperature and Hamiltonian of 
the system, respectively. 

Below, taking into account that the expression for the 
Green's function in the paraphase is independent of spin 
polarization, we omit spin indices. An important feature 
of the introduced functions is that I^(k, iujn) decomposes 
into the product of the propagator part and the strength 
operator- 
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(5) 



Thus, it is easy to obtain the modified Dyson equation 
for G(k, iun) 



ot ot 




(6) 

The bold fine in the equation corresponds to the total 
propagator G(k, iujn) and the triangle with symbol P 
denotes strength operator P(k, zco'n). The circle with in- 
scribed symbol corresponds to the Larkin-irreducible 
self-energy operator i;^(k, zcjn) The fine line with 
the light (dark) arrow denotes the seed Green's function 
for a Hubbard fermion that corresponds to the analytical 
expression 



Go{iuJm) 



1 



(7) 



The wavy lines with the light and dark arrows denote the 
Fourier image of hopping integral tk- The total propa- 
gator G(k, zcjn) relates to the strength and self-energy 
operators a; 

G(k,icjn) 
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i - t\^P{k, iujn) - S^(k, iujn) 



,(8) 



where ^ = e — ji. Making the analytical continuation 
iuJn 00 ^ i5 and introducing the real and imaginary 
parts of the strength and self-energy operators 

P(k,zcJn) ^i^(k,cj + i(5) = Pi(k,cj)+iP2(k,a;), (9) 



we arrive at 



P>(k,cj + id) = 



Pi(k,cj) + iP2(k,a;) 



^10) 



j-e-Ei(k,cj) + i((5-S2(k,a;); 

In this expression, 

Si(k,cj) = tkPi(k,cj) + (k,cj), 

S2(k,cj) = tkP2(k,cj) + S^(k,cj), (11) 

the real and imaginary parts of the Dyson-irreducible 
self-energy operator, respectively. 

Using the representation for the retarded Green's func- 
tion ([To]), we find the spectral intensity 

A{k,uj) -- Im D{k,uj) 

77 

1 [ [c.-g-Sf(k,c^)] •P2(k,u;) 
TT \ [w - e - Ei(k,w)]' + [S- E2(k,w)]' 



[5-E^{k,uj)]-Pi{k,uj) 
[w-C-Si(k,w)]' + [5-S2(k,w)]2 



(12) 



This formula establishes the interrelation between the 
spectral intensity ^(k,^^) and the strength P(k,co') and 
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self-energy i;^(k,a;) operators. Note that the denomi- 
nator of Expression ([12]) includes the Dyson self-energy 
operator, whereas the numerator includes only its Larkin- 
irreducible part. This resulted from mutual reduction 
of the numerator terms that are the product Pi(k,co') • 
P2(k, u). Our calculations demonstrate that the presence 
of non-zero imaginary part P2(k, uj) of the strength oper- 
ator in the numerator of expression (p!2]) leads to spectral 
intensity dependency of quasi-momentum on the Fermi 
surface. The pseudogap state of the normal phase of 
strongly correlated electron systems is attributed to this 
dependence. In the one- loop approximation^^^i^, the 
correction to P(k, iuJrn) caused by the interactions of the 
t — J model is determined by the four graphs 

w 



ti 



ot of 




ot 



ot 



(13) 



The contribution to the self-energy operator is deter- 
mined by the two plots 

of of of Q\ of of (14) 

In diagrams (13)-(14), the wavy line with the arrow de- 
notes hopping integral tq in the momentum representa- 
tion. The end of this line with the dark arrow forms 
the diagram fragment induced by the operator X^^ . The 
wavy lines without arrows denote the exchange integrals 
Jq. The longitudinal interaction Jfrn^^^^m shown 
by the wavy line with the two large circles. The end 
with the light circle corresponds to the diagram frag- 
ment in which the operator participated in pair- 
ing. The shaded circle corresponds to the operator X^^ . 
The transverse interaction JfmXJ^ X^ is denoted by 
the wavy line. At the ends of this line, sequence of two 
opposite values of the spin moment projection is shown. 
This sequence unambiguously points out that of the two 
operators describing the transverse interaction, pairing 
with which induced this diagram fragment. The dashed 
line corresponds to the Fourier image D_\_{q^,iuJi) of the 
quasi-spin transverse Green's function 

-(T.X^^(r)Xf (rO) = 



T 



— "^expii [q(R/ - Kg) -ujs{r -t)]} 



xP)^(q, zcjg). Us = 2s7rT, 



(15) 



and the shaded oval corresponds to the Fourier image of 



the Green's function 

- (t.A (X°°(r) + XjHr)) A (x™(r') + (r')) ) • 

In this expression, we used the notation 

A{Af{T)) = Af{T) - {Afir)). 

It can be easily seen that the function introduced for the 
oval can be expressed via the spin longitudinal Green's 
function and the charge Green's function. These two 
functions and their Fourier images are related as 



T 
N 



"^expii [q{Kf -Kg) -ujs{r -r')]} 



xL)||(q,zcjs), = 257rT, 



{TrAhf{T)Ahg{T')) 



(16) 



— exp{i [q(R/ - R^) - cj,(r - r')]} 



where 



(C(q, zcjs). Us = 2s7rT, 



An/ = X^^ + X^^ - {X^^ + X^^) , 



(17) 



= {X^^-X^^)/2, 

Associating analytical expressions to plots (p!3]) - (p!4|) . we 
obtain the strength and self-energy operators in the ex- 
plicit form 



P(k, iujm) = + ]^ X] (^q + -^k-q)G (q , iui) 
xx(q-k,zcj^ -ium), 

S'^(k) = - ^ ^ (tq + Jk-q)G (q, iui) , (19) 



(18) 



where Cn = 1 — n/2 is the Hubbard renormalization. 

These expressions and the faithful representation of 
the Green's function G (q ^iui) show that to solve the 
equation determining the strength operator, one should 
know the spin-charge susceptibility 

X(q, ium) = X5F(q, i^m) + XCF(q, ^^m), (20) 

that determines the contribution of the fluctuation pro- 
cesses. For convenience, this expression contains the dy- 
namic spin susceptibility 

X5F(q,^^m) = -D±{ci,iuJm) - D\\{ci,iuJm) 

= -3D\\{ci,iujm), (21) 
and the dynamic charge susceptibility 

XcF(q, i^m) = ^G(q, ium)' (22) 
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Expression (|2T]) was written taking into account the 
equality D±{q_^iujm) = 2D|| (q, zcj^), since without mag- 
netic field and the long-range magnetic order, the Hamil- 
tonian of the system is invariant relative to the transfor- 
mation of the SU{2) group^^'^^. 

Below, taking into account that the energy of charge 
excitations is relatively large, we limit our consideration 
to the contributions related to SFs. In this case, the 
nonlinear integral equation for the one-loop correction to 
the strength operator ^P(k, iuJrn) = ^(k, ioJrn) — Cn is 



N ^ ibJi 

Cl,lUJl 



^q-iq<^P(q,iw™)-S^(k)' 



(23) 



where = £ 
hopping integral is 



fi and the Fourier image of the 



= ^Iq + ^q, ^Iq = 2t{cOS Qj, + COS^^), 

i?q = At' COS Qx COS Qy + 2t'\cos2qx + cos2g'^). (24) 



III. MAGNETIC SUSCEPTIBILITY 

Since the kernel of the integral equation is determined 
by the dynamic magnetic susceptibility, let us briefiy an- 
alyze this function. Susceptibility xsri^^i^m) for the 
Hubbard model was calculated first in^^. Later, during 
the intense studies of high-temperature superconductiv- 
ity, xsri^^i^m) was calculated rpj^^ results of 
these studies show that the function Xsf{^.^ ^^m) rapidly 
drops with increasing Matsubara frequency. Therefore, 
the main contribution to the integral equation is collected 
by summation over ui close in value to cj^. Hence, we 
may assume that 



XSF (q, i^m) = xi^hsF ■ 



mO- 



(25) 



The expression for x(q) is the spin susceptibility at zero 
Matsubara frequency usf ^ ^sf/T^ where VtsF is the 
value of the Matsubara frequency starting from which the 
susceptibility rapidly drops. The order of magnitude of 
this frequency is determined by the characteristic values 
of excitation energies in a spin subsystem VtsF ^ 0.01 

It is important for further consideration that, in the 
weak doping region, the x(q) dependence in the t — J 
model is characterized by a sharp peak in the vicinity of 
the antiferromagnetic instability point Q = (7r,7r). Re- 
sults of the numerical calculations of the x(q) dependence 
using the technique from^^ are shown by the dashed line 
in Fig. [TJ They are in good agreement with the experi- 
mental data^^. 

To accelerate the numerical calculation in solving inte- 
gral equation (|34|) . we used the model susceptibility^^^^^: 



x(q) 



xo(0 



l + ^'(l + 7lq)' 



(26) 



where 



Xo(C) = 



3n 



4w.C(0 



, =0.55|t|, ^ = 2.7, (27) 

(cOS^a; + COSQy) 



e(l+7lq)^ 



Validity of this approximation follows from comparison 
of the dashed and solid lines in Fig. [H 




FIG. 1. (Color online) Quasimomentum dependences of spin 
susceptibilities: calculation by the t — J model according tc^^ 
(dashed line) and calculation by model susceptibility ([26]) 
(solid line). The chosen direction for rounding the Brillouin 
zone is r(0, 0) M(7r, tt) X{ti, 0) r(0, 0). 



With allowance for the above assumptions, we obtain, 
in the first Born approximation, that 



^P(k, iujr) 



1 ftq + -^k-q)x(q - ^)^SF 1^28) 



N 



^(^) = -^E(^q + ^k-q)nF(^q), 



(29) 



where np^x) = (exp(^^) + l)~-^ is the Fermi-Dirac func- 
tion. Then, expression (|5]) for the single-particle Green's 
function acquires the form 



Cn + 5P{\<i,iuJm) 



■a-M^(k,icj^)-S(k) 



.(30) 



The obtained system should be added with the equation 
for chemical potential fi 



n 



N 



E 



(31) 



The spectral intensity of the system can be calculated 
after analytical continuation using expression ([12]). 
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IV. PSEUDOGAP BEHAVIOR OF THE 
SPECTRAL INTENSITY OF HUBBARD 
FERMIONS 

Fig. [2]presents spectral intensities A(k, u) of the model 
under consideration calculated for the electron density 
n = 0.95 along the principle directions of the Brillouin 
zone (left plots) and on the Fermi contour for a quar- 
ter of the Brilouin zone (right plots). Here, the value 
of the spectral intensity is reflected by brightness of the 
energy spectrum lines (the lighter the line, the larger is 
the value ACk^uj)). The energy parameters of the model 
was measured in units \t\ 

T = 0.01, t = -1, = -0.65, t'' = -0.5, J = 0.25 

and chosen such that the Fermi surface would have the 
form of pockets, in accordance with the experimental 
data on magnetic oscillation^. The upper panel shows 
A(k, uj) corresponding to the Hubbard-I approximation. 
In this case, the value of A(k, uj) remains invariable at the 
change in the quasimomentum along the energy spectrum 
and at the Fermi surface. The middle panel demonstrates 
the results of the calculation of A(k,a;) with regard to 
SFs. Comparison with the upper panel shows that the 
allowance for SFs results in the qualitative difference, 
specifically, the occurrence of considerable A(k, uj) mod- 
ulation both at the spectrum line and at the Fermi level. 
It can be seen that the value of ACk^uj) decreases most 
in the wide energy region near the chemical potential. 

Note an important feature related to the effect of the 
sign of t on modulation of the spectral intensity. This 
feature is illustrated in the lower panel of Fig. [21 where 
the spectral intensity calculated for positive t is shown. 
In the calculation, all the rest parameters of the system 
remained invariable. Comparison of the middle and lower 
panels of the figure shows that at negative t the spectral 
intensity at the Fermi contour in maximum in the vicinity 
of the point (7r,7r), whereas at t > the A(k,a;) maxi- 
mum at the Fermi contour is located on the opposite side 
of the pocket, i.e., at the contour fragment that is closer 
to the point (0, 0). Note that this case corresponds to the 
results of the ARPES experiments (see, for example,^^). 

The developed modification of the spectral intensity 
by the expense of the fluctuation processes qualitatively 
changes the density of electron states 

^?H = ^E^(k,a;). (32) 

k 

Fig. [3] illustrates the densities of states calculated in 
the Hubbard-I approximation (dashed line) and with al- 
lowance for SFs (solid line). Comparing the two curves, 
one can see that the allowance for the fluctuation pro- 
cesses results in the drastic reduction of the density of 
states in the vicinity of the chemical potential. Thus, the 
results of our analysis show the formation of the pseudo- 
gap state in the considered system of Hubbard fermions. 




r M X r r x 




r M X r r 



FIG. 2. (Color online) Lines of the Fermi excitations and 
Fermi surfaces with allowance for the spectral intensity. The 
upper panel corresponds to the Hubbard-I approximation; the 
middle and lower panels, to allowance for SFs at different signs 
of hopping parameter t. 



V. SPECTRUM OF HUBBARD FERMION 
EXCITATIONS AND THE PSEUDOGAP 
BEHAVIOR IN THE LIMIT U ^ oo 

In this Section, we thoroughly analyze the physical na- 
ture of modulation of the spectral intensity in the system 
of Hubbard fermions. To elucidate the key points, we 
consider a simplified problem allowing us to obtain an- 
alytical expressions for the Fermi energy spectrum and 
spectral intensity but, at the same time, preserving the 
fundamental features of Hubbard fermions. It is possible 
with the use of the Hubbard model (j2j) in the regime of 
extremely strong electron correlations {U oo) when 
the Hamiltonian contains only the part corresponding to 
the operator of kinetic energy in the atomic representa- 
tion {t model) 

H = Ht. (33) 

Note that, since the model is simplified down to the limit, 
the below results on the properties of the Fermi excita- 
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FIG. 3. (Color online) Density of electron states calculated 
in the Hubbard-I approximation (dashed line) and with al- 
lowance for the spin- fluctuation processes (solid line). The 
vertical line reflects the position of the chemical potential. 



tion spectrum cannot be applied to the description of the 
experimental ARPES data and are used only to explain 
clearly the physical nature of the modulation of the spec- 
tral intensity in the ensemble of Hubbard fermions. This 
approach is analogous, in a sense, to the way in which 
the idealized Kronig-Penney model^^ illustrates the na- 
ture of the occurrence of energy bands in a crystal at the 
analysis of electron motion in a periodic potential. 

To study the spectral properties of Hubbard fermions 
in the t model, we use, as earlier, the method developed 
in Section [III In our case, the one-loop correction for 
^P(k, icj^) is determined only by the two upper plots 
from (p!3|) . The corresponding nonlinear integral equation 
for the strength operator will have the form 



SP{k,iuJrn) = T7 y] — 

N ^ luji 



^qX(q-k, iuji 



■ - tq (5P(q, iuJm) 



.(34) 



Choosing the quasimomentum dependence for the spin 
susceptibility, we again take into account the experimen- 
tal fact of the presence of a sharp peak of this function 
in the vicinity of the antiferromagnetic instability point 
Q = (7r,7r)^^. However, unlike the previous case, now we 
obtain simplified analytical expressions considering the 
limiting situation where the peak of magnetic suscepti- 
bility is delta-like. Then, Xspisi-,^) is presented as 



XSF (q,0) = xo • (^(q- Q). 



(35) 



To determine the value of xo, we use, as earlier, model 
susceptibility (|27|) . Then, for the first Born approxima- 
tion, we have 



(5P(k, iujm) 



iuJm — ^k+Q 



(36) 



and expression ([5]) for the single-particle Green's function 



acquires the two-pole structure 

Cn + (5P(k,icJm) 



P>(k, iujm) 



i^m - - tk(^P(k, iuJm) 
Cn{i^m - &+q) + ^5FXoWq 



(37) 



{lUJmY + ^'k^k+Q - ^^m(^k + Ck+q) " ^^SFXO^k^k+Q * 

Hence, the allowance for SFs forms two branches of the 
Fermi excitation spectrum 



^SFXoRi' 



The equation for chemical potential /i is 

lj2{nFiE+)A++npiE-)A, 



n 
2 



(38) 



(39) 



where the functions 



1 



-(Cn ± Ak), 
{C^ — ^SFXo)tl\i + ^SFXoR^k 

V{Cl-^SFXo)ti^^^SFXoIii 



(40) 



determine partial contributions of each branch of the 
spectrum to the total spectral intensity 

A(k,cj) = A+(k,a;) + A-(k,a;), 

As will be shown below, the change in the values of the 
functions at the change in k (for example, upon mo- 
tion along the isoenergetic line) leads to the formation of 
significant modulation of A(k, u) as soon as SFs become 
strong. 

The occurrence of A(k, uj) modulation is illustrated in 
Fig. m The upper plot on the left presents the energy 
spectrum development in the limit of the vanishingly 
small power of SFs {0.sf'Xo 0) for the electron density 
n = 0.95. As earlier, the values of the energy parameters 
of the model in the units of \t\ 

T = 0.01, t = -1, = -0.65, t'' = -0.4, (41) 

were chosen such that the Fermi surface would have the 
form of the pockets. The dashed hue corresponds to the 
lower branch of the spectrum; the solid line, to the 
upper branch . Despite there are two solutions of the 
dispersion equation, in reality the only branch is revealed. 
This is due to the fact that at ^sfXo the partial 
contributions 



= ^^n{l±sign(tik)}. 



(42) 



sharply grow from zero to the end value or sharply drop 
from the end value to zero at the variation in the quasi- 
momentum (dashed lines in Fig. [5]). Therefore, at each 
nonzero point of the Brillouin zone having the same 
value, only one of the partial spectral intensities A^ (k, uj) 
may exist. The total A(k,cj) spectral intensity (the val- 
ues of these functions are reflected in the plots by thick- 
ness of the corresponding lines) remains invariable at 
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FIG. 4. (Color online) Spectrum of the Fermi excitations 
along the principal directions of the Brillouin zone and modi- 
fication of the spectral intensities with allowance for the spin- 
fluctuation processes. 




FIG. 5. (Color online) Dependences of partial contributions to 
the spectral intensity on the quasimomentum with allowance 
for SFs (solid lines). For comparison, the dashed lines show 
the dependences of these functions in the Hubbard-I approx- 
imation. 
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the variation in the quasimomentum along the energy 
spectrum. In our simple case, this corresponds to the 
Hubbard-I approximation. The considered features are 
illustrated by the left plots of the second, third, and 
fourth panels in Fig. HJ 

The situation becomes qualitatively different when SFs 
switch (right plots in Fig. S]). The upper right plot 
presents the spectra calculated at ^sfX^ — 0.16. As 
before, the solid line shows the spectrum; the dashed 
line, the spectrum E^. Switching of the interaction be- 
tween the electron subsystem with the subsystem of spin 
degrees of freedom via the kinematic mechanism leads to 
repulsion of the spectrum branches at the points where 
they were osculating without interaction. 

The more significant effect of SFs is related to con- 



FIG. 6. (Color online) Effect of the sign of hopping parameter 
t on modulation of the spectral intensity. 



siderable modification of the spectral weights Aj, which 
manifests itself in strong renormalization of the values of 
these functions (Fig. [5|). As a result, the partial spectral 
intensity A+(k, cj) acquires a finite value over the entire 
curve of the energy spectrum E^ (similarly, A~ (k, uj) is 
finite along the entire E^ curve) and becomes strongly 
modulated. It can be clearly seen from comparison of 
the left and right plots in the second and third panels 
in Fig. m The total spectral intensity A(k, uj) acquires 
the structure reflecting the presence of two branches ex- 
hibiting strong modulation. In the right plot of the lower 
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panel in Fig. IH the degree of darkening of the parts in 
the vicinity of the points (k, uj) corresponds to the value 
of A(k,a;) in these points. It can be seen that the al- 
lowance for SFs has led to induction of the shaded zone 
and redistribution of the spectral intensity between the 
basic and shaded zones. 

The feature related to the change in the sign of param- 
eter t of hopping to the first coordination sphere (Figl2j) 
is characteristic of this model. Fig. [6] shows the spectral 
intensity calculated along the principle directions of the 
Brillouin zone (left plots) and at the Fermi contour for 
a quarter of the Brillouin zone (right plots) at different 
signs of t. It can be seen that, at t > 0, ^(k,^^) is max- 
imum in the part of the Fermi contour that is closer to 
the point (0, 0), which is consistent with the experimental 
data. 



VI. SUMMARY 

To sum up, we formulate the principles of the kine- 
matic formation of modulation of the spectral intensity 
A(k^uj). Of fundamental importance is the use of the 
faithful representation for a single-fermion Matsubara 
Green's function DCk^iujm)- For the kinematic mecha- 
nism, D(k, iujm) is expressed as the product of the prop- 
agator part and strength operator P(k, zcj^n)- The pres- 
ence of the strength operator in the numerator of the 
Green's function and its dependence on the Matsubara 
frequency and quasimomentum lead to the fact that the 
isoenergetic lines in the quasimomentum space become 
different from the lines where the strength operator has 
a constant value. This spacing is one of the causes of 
modulation of the spectral intensity A{]<.,uj). It is impor- 
tant that the integral of hopping to the first coordination 



sphere determines the Fermi contour part where A(k, uj) 
considerably decreases. 

The specific cause of A(k, uj) modulation in the frame- 
work of the kinematic interaction is that SFs lead to the 
formation of a shaded zone, which represents the initial 
zone shifted in the quasimomentum space by the vector 
(7r,7r). As a result, the total pattern of the spectrum is 
formed by coherent hybridization of these two zones, be- 
tween which the spectral intensity is redistributed and, 
consequently, the density of states in the vicinity of the 
chemical potential drops. 

To demonstrate the kinematic formation of the pseu- 
dogap behavior as brightly as possible, we considered the 
Hubbard model in the limit of strong correlations when 
the dynamics of Hubbard fermions is governed by the 
kinematic interaction. Obviously, the discovered mecha- 
nism of the pseudogap phase formation is universal and 
relevant for other models of strongly correlated systems. 
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